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The static friction between crystalline surfaces separated 
by a molecularly thin layer of adsorbed molecules is calcu- 
lated using molecular dynamics simulations. These molecules 
naturally lead to a finite static friction that is consistent with 
macroscopic friction laws. Crystalline alignment, sliding di- 
rection, and the number of adsorbed molecules are not con- 
trolled in most experiments and are shown to have little ef- 
fect on the friction. Temperature, molecular geometry and 
interaction potentials can have larger effects on friction. The 
observed trends in friction can be understood in terms of a 
simple hard sphere model. 



I. INTRODUCTION 

The last decade has seen great advances in techniques 
for measuring friction in contacts whose geometry And /or 
chemistry are controlled with atomic precision.lJT3 At 
the same time, increases in computational power hayje 
allowed increasingly realistic simulations of frictionEm 
These studies reveal unexpected behavior that raises 
questions about the molecular origins of static friction 
and the "laws" of friction that describe most macroscopic 
objects. 

Friction is the lateral force needed to slide one object 
over another. The force needed to initiate sliding is called 
the static friction F s . Its existence implies that the sur- 
faces have locked into a local free energy minimum, and 
F s represents the lateral force needed to displace them 
out of this minimum. The force needed to maintain a 
constant sliding velocity v is called the kinetic friction 
Fk{v), and it represents the force required to replace en- 
ergy dissipated during sliding. 

Macroscopic objects almost always exhibit a finite 
static friction and a kinetic friction that is slightly smaller 
at low velocities. One puzzling result from many molec- 
ular scale theories of friction between bare surfaces is 
that the static friction almost always vanishes, and is 
not closely related to the kinetic friction.Lrcl This indi- 
cates that some important feature is missing from these 
model systems, that must be included to make contact 
with macroscopic experiments. Of course real surfaces 
are generally not bare, but are coated with a layer of 
adsorbed molecules, as well as dust, dirt and other de- 
bris. In this paper we explore the influence of a layer 
of molecules between two surfaces on friction forces, and 
show that including such layers naturally leads to behav- 
ior that is consistent with macroscopic measurements. 

Over three hundred years ago Amontons' published 
two "laws" of macroscopic friction that are still taught 



and used today.0 These state that the friction is pro- 
portional to the normal load L pushing two surfaces to- 
gether, and independent of the apparent geamefarical area 
of these surfaces A app . Subsequent workBaOll3 is con- 
sistent with a more_general friction formula that is based 
on the observatiorJlZTLj that the actual area of molecular 
contact between two surfaces A Tea \ is much smaller than 
the nominal surface area ^4 a pp- The friction is assumed 
to be given by ^4 roa i times a local shear stress r. If r rises 
linearly with the local contact pressure P, 



T = T + CtP , 

then since P = L/A lea \ 

F = rA Ieal = T A rcal + ah 



(1) 
(2) 



This expression agrees with Amontons' laws if A rea i is 
proportional to L or if To is sufficiently small compared 
to aP. The former condition applies if the load is 
high enough to produce plastic deformation of thc_siii'- 
faces so that L/A Tea \ remains equal to the hardnessO It 
also holds for non-adhesive contacts betwEep_.ideal elas- 
tic solids with random surface roughness. cj£j However, 
any adhesive bonding leads to friction in the limit of zero 
load and violates Amontons' laws. For a piece of adhesive 
tape the first term in Equation ^ dominates at low loads, 
and in many cases friction can be observed at negative 
loads. 

Early attempts to explain Amontons' laws and the 
existence of static friction were based on the idea that 
peaks cui .one surface interlock with valleys on the other 
surface £30 In order to slide, the top surface must then 
be lifted up a ramp formed by the typical slope tan ip of 
the bottom surface. If there is no microscopic friction 
between the surfaces, then the minimum lateral force 
to initiate sliding is F s = Ltanip. This result satis- 
fies Amontons' laws with a constant coefficient of fric- 
tion )j, s = F s /L = tan^> that can span all possible val- 
ues. However, this geometrical model for friction can 
not explain many experimental observations. For exam- 
ple, coating a single monolayer of surfactant on a surface 
does not change its slope, but can reduce the friction by 
more than an order of magnitude. It is also well known 
that making surfaces too smooth actually leads to an in- 
crease in friction.E3E3 A practical illustration of this is 
that magnetic hard disks are purposely roughened to re- 
duce friction. 

The fundamental problem with the above explanation 
for the origin of static friction is that there is no reason for 
the peaks from different surfaces to be correlated. In gen- 
eral one expects that at any instant in time some peaks 
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will be moving up a ramp and an equal number of peaks 
will be moving down. The net lateral force will average to 
zero and there will be no static friction. A similar prob- 
lem arises if one imagines that the interlocking peaks are 
individual atoms as in the "cobblestone" model discussed 
by Israelachvili and coworkers. Ejo If atoms are from two 
identical, aligned surfaces all of them will go up ramps 
simultaneously yielding Amontons' laws. However, fric- 
tion measurements are usually made between misaligned 
surfaces or surfaces with different lattice constants. In 
this case the force from atoms moving up ramps should 
cancel the force from atoms moving downward, yielding 
a vanishing static friction. ■— ■ ■— . 

A nuinber—. of detailed analyticcilEa and 
simulatiorJ3n3i2fl studies have concluded that F s should 
generally vanish between bare crystalline or disordered 
surfaces. The exception is the case of commensurate sur- 
faces which share a common periodicity in their plane 
of contact. In this case the free energy varies with the 
phase difference between the common—Eaurier compo- 
nents, and simulations of such surfacesEjtj find a static 
friction that rises linearly with pressure as in Eq. [I]. 
However the probability that two contacting surfaces are 
commensurate is infinitesimally small. Even identical 
surfacesi-are only commensurate when they are perfectly 
aligned.Ea Any misorientation (Figure 0) causes the sur- 
faces to become incommensurate, i.e. have no common 
period. The free energy is thcpjpdepcndent of transla- 
tions and F s is identically zeroE3E3'E3 If the interactions 
between the two surfaces are strong enough compared to 
interactions within the surfaces, they may deform into a 
commensurate structure.Cj'EErtj However, this criterion 
is unlikely to be met for most of the surfaces around us, 
which have been chemically passivated by exposure to 
air. Indeed calculations for clean surfaces in vacuum find 
zero static -friction between different facets of the same 
metal£30'E3E3 Moreover, contacts that ara-kpawH to be 
in the elastic limit exhibit static friction.uoEjlIa Edge 
effects can lead to a finite static friction, but are not 
significant for the micron size-contacts typically found 
between macroscopic objects.L3 Chemical disorder and 
surface roughness, axe also unable to produce observed 
friction forces EJHcj 

There are relatively few direct experimental measure- 
ments of the friction between bare crystalline surfaces. 
However, stedies otadsorbed monolayers, c2lc3 smalLcrys- 
talline tips,c2l mica,E3 and molybdenum disulphideca sup- 
port the conclusion that friction vanishes in this limit. A 
finite friction has been observed between clean metallic 
surfaces, EZrc3 but the surfaces were rough and the load 
appears to have been high enough to produce plastic de- 
formation. Further work in this area would be of great 
value in testing the. abijOffi-theoretical models. 

In recent papersOuO'ufl we have suggested that the 
so-called "third bodies!! that are present between most 
contacting surfacesou can provide a general mecha- 
nism for static friction. These third bodies may take the 
form of wear debris, dust, or small hydrocarbon or water 



molecules that are adsorbed from the air. The interac- 
tions between such bodies are generally weaker than the 
interactions within the bounding solids. This frees them 
to rearrange at the interface to lower their free energy 
and lock the surfaces together. The layer of third bod- 
ies creates an immense number of metastable states, like 
that in a glass or granular medium. The third bodies can 
always fall into one of these states that is in registry with 
both bounding walls and thus produce static friction.Eil 

In this paper we report extensive studies of the static 
friction produced by layers of spherical or short chain 
molecules between crystalline surfaces. We begin by de- 
scribing different methods of measuring the static fric- 
tion and establishing that it has a well-defined thermody- 
namic limit. Then the effect of temperature, wall geom- 
etry, interaction potentials, chain length and areal den- 
sity on the static friction are explored. In each case the 
static friction obeys Equation [I] over the experimentally 
relevant pressure range. Factors that are not controlled 
in typical macroscopic experiments have little influence 
on To and even less effect on a, which dominates the fric- 
tion at high loads. Such factors include wall geometry, 
sliding direction, and the length and density of adsorbed 
molecules. Increasing the temperature lowers To, but has 
relatively little effect on a. The value of a is strongly 
dependent on the interaction potential, particularly the 
effective hard sphere size of the molecules compared to 
that of wall atoms. All of our results can be understood 
in terms of a simple geometrical explanation for the ori- 
gin of Equation [y. 

The following section describes the details of our simu- 
lations and the algorithms used to determine F s . Section 
III determines the effect of each parameter in our simu- 
lations on the static friction and describes a simple geo- 
metrical picture that explains all of the observed trends. 
The final section summarizes the results. 



II. SIMULATION METHOD 

A. Potentials and geometry 

We use a bead spring model@ that allows us to ex- 
plore the behavior of simple spherical molecules or short 
chains between bounding solids. Each molecule contains 
n spherical monomers. All monomers interact with each 
other through a truncated Lennard-Jones potential: 



V(r) = 



4e[((j/r) 12 




{a/rf - (a/r c ) 12 + (<r/r c ) 6 ] ; r < r c 

;r > r c 

(3) 



The parameters a and e are chosen as the units of 
length and energy, respectively. Combined with the 
monomer mass m, they determine the unit of time 
£lj = (raa 1 1 1) 1 ! 2 . Typical values for hydrocarbons are: 
a ~ 0.5nm, e ~ 30meV and £lj ~ 3ps. We consider 
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the case r c = 2 1 / 6 er to simulate a purely repulsive poten- 
tial, and r c = 2.2a is used to study the effect of adhesive 
forces. 

Monomers are bound into chains by an additional 
strongly attractive potential between nearest neighbors 
on the same molecule: 



V 



<:h m = f -ikRoMl - {r/Ro) 2 }, r < R 
( ' I oo r> R 



(4) 



where i?o = 1.5c and k = 30ecr~ 2 . These parameters are 
chosen to_Drevent chain crossing as described in previ- 
ous work.E3 We consider the case of spherical molecules 
(ri = 1), and of short chains with n = 3 and 6. Based aa 
mappings of the bead-spring model to real molecules ,lj 
this corresponds to alkane chains with up to roughly 16 
carbons. Substantially longer molecules are unlikely to 
have sufficient volatility to contribute to the airborn con- 
tamination that is present on any surface. Polymers are 
often placed at interfaces intentionally to act as lubri- 
cants. We have not addressed the behavior of these much 
longer molecules. 
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Molecules are confined between two walls formed by 
the (111) surfaces of fee crystals directed normal to the 
z-axis. Wall atoms are connected to their lattice sites 
with springs of stiffness k s in ocder to model elastic de- 
formation in the simplest wayS We consider the com- 
pletely rigid case, k s = oo, and k s — 840e/<7 2 or 210e/er 2 . 
Both of the latter values correspond to relatively com- 
pliant solids. For example, at T = 0.8e/fcs they lead 
to root mean squared (rms) displacements about lat- 
tice sites that are about 4.5 and 9%, respectively, of the 
nearest neighbor separation (1.2<r) used in most of our 
simulations £3 

Wall atoms and fluid monomers also interact with a 
Lennard-Jones potential but with different characteris- 
tic energy and length scales er w f and e w f, respectively. 
These are varied to determine the effect of molecular size 
and chemistry on static friction. Direct interactions be- 
tween atoms on different walls are not included in the 
simulations. They vanish identically for the short cutoff 
r c = 2 1 / 6 er used in most of our simulations, and have 
a very small effect on calculated quantities at larger r c . 
Tests of the effect of wall-wall interactions are discussed 



briefly in Section III F 



Periodic boundary conditions are imposed in the plane 
of the walls. These periodic boundary conditions prevent 
us from considering truly incommensurate systems. How- 
ever, the effect of commensurability r apidly-d ffirea ses as 
the length of the common period increases EaSE3 The 
degree of commensurability between the walls is varied 
in two ways. The first is to rotate the top wall by an an- 
gle 9 about the z-axis (Fig. |A-C). Only the range of 9 
from to 30° produces inequivalent results. The second 
is to retain the alignment of the two lattices (9 = 0), but 
to change the ratio of the lattice constants of the bottom 
and top walls from unity to 12/13 (Fig. |l| D). 
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FIG. 1. Projections into the x — y plane of atoms from 
the inner surfaces of the bottom (solid circles) and top (open 
circles) walls. In A through C the two walls have the same 
structure and lattice constant, but the top wall has been ro- 
tated by 0°, 8.2°, or 90°, respectively. In D the walls are 
aligned, but the lattice constant of the top wall has been re- 
duced by 12/13. The atoms can only achieve perfect registry 
in the commensurate case A. The simulation cell was usually 
at least four times the area shown here. 



TABLE I. Values of the orientation angle 6 and the corre- 
sponding two dimensional basis vectors and number of atoms 
per surface layer in the top wall. Both surfaces have hexagonal 
symmetry. The final column gives the percentage difference 
between the nearest-neighbor spacings in the two walls. 



6 


bi/d 


^ Id 


No. of atoms 


(%) 


0° 


i 


2 l + 2 J 


2304 





2.02° 


588 7* , 12V3 7/ 
601 ' 601 J 


276 7* I 300 V^T* 
601 ' 601 J 


2404 


-2.1 


4.127° 


192 7/ , 8^3 7" 

193 ' ~ 193 J 


84 7* 1 100^3 7; 
193 ' 193 J 


2316 


-0.3 


6.3° 


564 7* , 36V3 7/ 
559 559 J 


■22*7 1 300 x/3 ? 
559 T* 599 ■> 


2236 


1.5 


8.2° 


487, 

49 ' ~ 48 J 


187 , W37; 

49 1 t 4S ,y 


2352 


-1.03 


10.4° 


564 7 1 , 60v^7 
571 ' ~ 571 J 


192 7* 1 312^3 7; 
571 ' ~ 571 J 


2284 


0.44 


14.4° 


564 7 , 84^3 7/ 
589 ' ~ 589 J 


156 7 1 324 v 7 ! 7 
5S9 1 589 J 


2356 


-1.1 


19.1° 


20 7 , 4^3 7 

21 21 J 


21 1 + 21 -< 


2268 


0.79 


25.2° 


33 7 , 9%/3 7 
.'17 1 ^ ZJ_J 


3 7 , 21V37 

37 4 + 37 J 


2368 


-1.4 


30°/90° 




4^3 7/ 
7 J 


2352 


-1.03 
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The simulation cell is a rectangle whose height along y 
is v3/2 times the length along x, and the [110] direction 
of the bottom wall is directed along the x axis. Special 
values of 9 are chosen that allow both top and bottom 
surfaces to retain perfect triangular symmetry, and to 
have nearly the same nearest-neighbor spacings. Values 
for the percentage difference between the spacing on the 
bottom, d, and on the top, d! , walls are given in Table 
The difference is typically less than 2%, and no trends 
were seen with the sign or magnitude of this difference. 

The system is thermostatted by coupling to a heat 
reservoir through Langpifin noise and damping terms in 
the equation of motionES Previous work shows that cou- 
pling along the direction of sliding produces a direct effect 
on the kinetic friction, while .coupling to the perpendicu- 
lar components does not.l3l23'L3 Thus only the perpendic- 
ular components are thermostatted, using the equation 
of motion: 



s 



mx± = Fj_ — mTx±_ + W± 



(5) 



Here F is the force from the interaction potentials, T 
is the friction constant that controls the rate of heat 
exchange with the reservoir, and W is the Gaussian- 
distributed random force acting on each monomer. The 
rms value of W is determined fpam T through the usual 
fluctuation-dissipation relation.E3 The equations of mo- 
tion are integrated insing a fifth-order, Gear predictor- 
corrector algorithm. In most cases the time step dt = 
0.005i L j and V = 0.4*7? . As discussed in Sec. |ffl A|, 



our measured static friction forces are insensitive to the 
choice of T. 

The areal density of monomers is specified as the ef- 
fective coverage on the separate surfaces before they are 
brought into contact. To avoid ambiguities in the defi- 
nition of the monolayer density, we define the coverage 
as the number of adsorbed monomers divided by the to- 
tal number of wall atoms in the surface layers of both 
walls. A coverage of 1/2 means that each wall had half 
as many adsorbed monomers as surface atoms before the 
walls were placed in contact. After contact the monomers 
from the two walls will produce about a monolayer be- 
tween the walls, because the monomers and wall atoms 
have roughly the same size in our simulations. 

To prepare initial configurations, we typically start 
with monomers in a crystalline state, or in the final state 
of a simulation for different conditions. The system is 
then heated to a temperature of 1.9e/ 7c_b at a constant 
pressure of 4e<7~ 3 until the configuration has randomized 
(typically 500iL,j). Finally, the temperature and pressure 
are ramped to the desired value over 50 £lj and allowed 
to equilibrate there for 450£lj ■ Results for the static fric- 
tion are not sensitive to changes in this procedure. 



B. Determining the static friction 

As noted in the introduction, static friction arises be- 
cause the system has managed to lock into a local free en- 



ergy minimum. The total energy needed to activate the 
system out of this free-energy minimum increases with 
increasing system size.Ea However, for any finite system, 
thermal fluctuations will eventually lead to activated dif- 
fusion in the absence of any lateral force. This diffusion 
has been studied previously for the same model system 
used here.E-3 As discussed by these authors, the static fric- 
tion is only well-defined if the limit of infinite system size 
is taken before the limit of long times. If the limits are 
taken in this order, then the static friction corresponds 
to the maximum derivative of the free epergy as the sys- 
tem is driven out of its local minimum EB If the order is 
reversed, the static friction is always zero. 

Note that even in the infinite-size limit there may be 
rate dependence in the measured static friction. This 
is because the effective free energy surface depends on 
the degree to which monomers can diffuse between the 
walls. Based on our recent studies of kinetic friction, we 
expect that monomer diffusion will lead to a weak log- 
arithmic dependence of static friction on measurement 
timeO This is analogous to the logarithmic rate depen- 
dence in the apparent yield stress of glassy systems.EJ 

Two methods are used to find the static friction. We 
call them "ramp" and "search" respectively. In both 
cases, simulations are done with the bottom wall fixed 
and the top wall under constant normal pressure P. A 
lateral or shear stress r is then applied to the top wall at 
an angle (f> relative to the x-axis. In the ramp algorithm, 
the shear stress is increased from zero at a constant rate 
until the wall begins to slide. The stress at which motion 
initiates is recorded, and the stress is then decreased at 
the same slow rate. The wall stops moving as the force 
decreases towards zero, and then begins to slide again 
in the other direction when the force is sufficiently nega- 
tive. The magnitude of this threshold force is recorded, 
and the stress is increased once more. This oscillatory 
process is repeated at least 15 times to get a statisti- 
cal sample of threshold forces. The stress is typically 
changed at a rate of 0.002e/eriLj- This rate must be slow 
enough that the system has time to repin as the mag- 
nitude of the force drops. It also sets the accuracy on 
determining the depinning force, because a finite time is 
needed before motion of the top wall can be detected. In 
most cases the force is increased until the wall reaches a 
relatively high sliding velocity (0.5ct/£lj) and the depin- 
ning point is identified with the end of the last previous 
interval of 20£lj during which the wall moved less than 
0.5ct. 

In the search algorithm, we first guess upper and lower 
bounds for the static friction. Then a lateral stress equal 
to the average of these, values is applied to determine 
whether sliding occurs.EJ The top wall is considered to be 
stuck if it moves by less than la in lOOiLj and to be slid- 
ing if it moves by more than 5c. If it moves an interme- 
diate distance, the motion is followed for another 100£lj, 
and the wall is considered to be sliding if it moves by more 
than it did during the first 100£lj- If the trial stress does 
not lead to sliding it becomes the new lower bound, and 



4 



if sliding does occur the trial stress becomes the new up- 
per bound. The system is then returned to an identical 
initial condition, and a stress midway between the two 
new bounds is applied. This process is repeated enough 
times (typically five) to obtain the threshold stress with 
a precision of better than ±0.02ec~ 3 . To get good statis- 
tics, different initial configurations are tried. These are 
generated by allowing sliding to occur for a distance of at 
least 5<T and then equilibrating the system at zero lateral 
stress for at least IOO^lj- 



III. RESULTS 

We first consider how technical details of the simula- 
tions may influence the results. These include the effects 
of thermostatting, system size, and the choice between 
ramp and search algorithms. Then the relevant physical 
parameter space is explored. First the variations with 
pressure, temperature and wall stiffness are discussed. 
Then the effects of parameters that are not controlled in 
experiments are considered. These include the wall ori- 
entation 9, sliding direction <j), coverage, and chain length 
n. Finally, the effect of the potential parameters cr w f , e w f , 
and r c is explored. 



A. Effect of thermostat 

As described in the previous section, the wall atoms 
and monomers are coupled to a thermostat in order to 
maintain a constant temperature. The rate of heat trans- 
fer r is chosen to be small enough that the thermostat 
has little effect on the dynamics over the time scale re- 
quired for the velocity autocorrelation function to decay 
to zero. 

The main function of the thermostat is to remove heat 
generated by sliding friction. Previous work shows that 
the thermostat has little effect on kinetic friction if it is 
only coupled to the velocity corajtanents that are perpen- 
dicular to the sliding direction.ErEH Our focus here is the 
static friction. This should be even less sensitive to the 
thermostat since no heat is generated until sliding starts, 
and this is at a force greater than the static friction. 

Fig. H compares results for the static friction as a func- 
tion of pressure P for V = 0.4 and l.QtjTj. The cov- 
erage was 1/2 and the other parameters for these sim- 
ulations were 9 = 90°, e w f = e, c w f = a, r c = 2 1 / 6 u, 
and T = 0.7e/fcs. Note that results for the two ther- 
mostatting rates are equal within our statistical error 
bars. Other tests show that equivalent results are ob- 
tained when only the wall atoms are coupled to the ther- 
mostat. 
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FIG. 2. Yield stress vs. pressure for two values of the pa- 
rameter F that controls the rate of heat exchange with the 
thermostat (Equation ||). We observed no influence of the 
thermostat on the yield stress. For this data, the coverage is 
1/2, 9 = 90°, <r wf = I.Oct, e wf = l.Oe, r c = 2 1 / 6 <r, d nn = 1.2a, 
T — 0.7e/fes and n—6. The search algorithm was used, and 
the simulation cell was 28.8a by 24.9<r. Error bars are indi- 
cated when they are bigger than the symbol size. 



B. Effect of system size 

The phenomenological expressions for friction de- 
scribed in the introduction assume that the static friction 
is given by a yield stress times the area of contact. The 
regions of direct molecular contact between macrescaruc 
solids typically have diameters of order microns .ETO'tS 
It is important to establish that yield stresses calculated 
at the much smaller scales of our simulations are rep- 
resentative of larger contacts. It is also interesting to 
consider how friction forces would change if the contact 
dimensions were reduced to nanometer scales. 

Figure |^ compares results from the ramp and search 
algorithms and shows how both scale with contact size. 
Values of the yield stress r s at P = 12 and 28ea~ 3 are 
plotted as a function of the inverse width of the simu- 
lation cell. The latter is expressed as A -1 / 2 where A is 
the area in the x-y plane. Given a ~ 0.5nm, the largest 
contact widths (smallest values of A -1 / 2 ) correspond to 
50nm and the smallest widths to about lOnm. 

We first consider how contact area affects the varia- 
tions in the stress needed to initiate sliding from differ- 
ent initial configurations. Points indicated by stars and 
squares indicate the maximum and mean forces needed to 
initiate sliding for P = 28ea~ 3 with the ramp algorithm. 
Note that the difference between the extreme and aver- 
age yield stresses decreases rapidly with system size and 
is about 10% for the largest systems. This indicates that 
the distribution of yield stresses is well-behaved and can 
be thought of as arising from an incoherent sum of contri- 
butions from different regions of the contact. The same 
trend is seen for the maximum and mean values from 
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the search algorithm, but the fluctuations are slightly 
smaller. 

In both algorithms the smallest system sizes tend to 
produce larger mean stresses. However*, as assumed in 
phenomenological theories of frictionJlJ'EjO the mean 
values of the yield stress for all pressures and algorithms 
go to well defined limits as the contact size diverges. As 
shown in Fig. |§|, finite size effects are larger for higher 
pressures. This may be because the energy barrier for 
monomers to move between metastable states increases 
with pressure, making them more likely to get stuck in 
high energy states. Finite-size effects also grow rapidly 
as the coverage decreases, since there are then fewer 
monomers at the same wall area. 

Note that yield stresses from the ramp algorithm tend 
to be consistently higher than those from the search al- 
gorithm. At P — 28e<7~ 3 the difference appears to reflect 
larger finite-size effects for the ramp algorithm. However 
the difference between the two algorithms is nearly in- 
dependent of system size for P = 12ea~ 3 . The reason 
for this discrepancy is that it takes time for the wall to 
accelerate once the static friction is exceeded. One can 
only be sure that the wall has started sliding when it has 
moved by a few a. With the ramp algorithm the force 
has increased to a higher value by this point. Attempts 
to extrapolate back to the onset of motion only partially 
corrected for this offset. Decreasing the ramp rate re- 
duces the error, but also lengthens the simulation time 
proportionately. The search algorithm reduces the effec- 
tive ramp rate to zero and efficiently converges on the 
yield stress. 
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FIG. 3. Mean yield stress as a function of the inverse square 
root of the area A of the periodic simulation cell for the 
ramp (squares) and search (circles) algorithms at P = I2ea~ 3 
(closed symbols) and P — 28ecr~ 3 (open symbols). The max- 
imum stress needed to initiate sliding from any initial con- 
dition at P = 28ecr -3 is indicated by asterisks for the ramp 
algorithm. Simulations were done with rigid walls at 1/2 cov- 
erage, 6 = 90°, (Twf = a, e w f — e, r c — 2 1 ^ 6 a, d nn = 1.2cr, 
T = 0.7e/fcs and n — 6. Statistical errors are indicated by 
the symbol size. 



For the above reasons the search algorithm is used for 
most of the following plots. Except where noted, the sec- 
ond to smallest system size of Figure [3] is used. The sur- 
face of each wall then has 576 atoms and the simulation 
cell is 28.8cr by 24.94a for the usual case of d nn = 1.2a. 
Larger system sizes are used mainly at low coverage in 
order to ensure that the deviation of the yield stress from 
the thermodynamic limit is comparable to statistical un- 
certainties. Larger systems also allow a greater number 
of rotation angles, 6, to be studied. 



C. Effect of temperature and pressure 

Figure |I] shows the yield stress vs. pressure for a 
coverage of 1/2 over a broad temperature range from 
ksT/e = 0.3 to 1.3. For reference, the triple point and 
critical temperature of Lennard-Jones monomfr.s with 
r c — ► oo are about 0.7 and 1.3 e/fce, respectively,E3E3 and 
the zero pressure glass transition temperature for—short 
Lennard-Jones chains is between 0.4 and 0.5e//cBO~l£3 
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FIG. 4. Variation of yield stress with pressure at the in- 
dicated temperatures. In each case, r s rises linearly with P. 
The size of typical statistical errors is indicated by the symbol 
size. Errors are slightly larger at the highest pressures. Simu- 
lations were done with the search algorithm for rigid walls, 1/2 
coverage, 9 = 90°, a w f = a, e wt = e, r c = 2 1//6 (j, dnn = 1.2a, 
and n=6. 



Note that t s increases linearly with pressure at each 
T, as assumed in Equation |l|. For representative val- 
ues of e=30mcV and a = 0.5nm, the unit of pressure 
eer~ 3 ^40MPa. Thus the highest pressures in Figure 
H correspond to about 1.5GPa. In those cases that we 
checked, linearity extended up to much larger pressures 
(at least lOOeer -3 ). However, these pressures were not 
routinely studied because a smaller time step is needed 
to integrate the equations of motion at larger pressures. 
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D. Effect of wall stiffness 



k B T/e 

FIG. 5. Values of a and to as a function of temperature 
from fits of Equation ^ to the data of Figure ^. Dashed lines 
show the linear fits to the temperature dependence given in 
Equation H. 



The plots of yield stress vs. pressure at different tem- 
peratures in Figure ^ are nearly parallel, but their inter- 
cepts shift rapidly with temperature. Results from fits 
to Equation [l are shown in Figure H The slope a drops 
by only 16% with increasing T, while the intercept To is 
roughly proportional to T. Linear fits as a function of 
temperature give: 



a 



= (-0.0057 ±0.0007)^ 
= (-0.214 ±0.013)-^%- 



+ (0.0338 ±0.0006) 
(0.015 ±0.011). 

(6) 



An explanation fo r the nearly constant value of a is pre- 
sented in Section III F . In the remainder of this subsec- 



tion we focus on the variation of To- 

At each T, the linear fit for yield stress vs. pressure 
crosses through zero at a threshold pressure Pt = —To/a. 
Above this pressure we find a well defined yield stress in 
the thermodynamic limit. Below P t , the adsorbed layer 
acts like a lubricating liquid, with a._&hear stress that 
vanishes as the shear rate decreases. e£| In this regime, 
the bounding walls are far enough apart that molecules 
diffuse freely and do not lock the two walls in a local en- 
ergy minimum. From the fit parameters in Fig. || we see 
that Pt rises roughly linearly with T. This is because the 
simulations used r c — 2 1 / 6 cr, which corresponds to the 
hard sphere limit. In this regime pressures are balanced 
only by entropic forces and thus must scale linearly with 
temp eratu re. Attractive interactions are considered in 
Sec. En? 



Most of the following simulations were done at T = 
0.7e/fcs- Given the results just described, we expect 
that the same trends would be observed at other tem- 
peratures. 



As described in Section II A, wall atoms were con- 
nected to lattice sites by springs with stiffness k s , or fixed 
rigidly to these sites (k s — > oo). Decreasing k s has two 
competing effects on the static friction. One is to make it 
easier for the two surfaces to deform into a structure with 
a common period. This tends to increase the yield stress. 
However, studies with bare walls found that unrealisti- 
cally small values of k s were needed to produce a finite 
static frictiono Decreasing fc s also decreases the yield 
stress by allowing wall atoms to be pushed out of the 
way. In our simulations the second effect always domi- 
nates, and r s decreases as k s decreases. This is consistent 
with a recent analytic theory for static frictionEfj 

Fig. |TJ shows results for systems C and D of Fig. |l] at 
ksT = 0.7e. In each case, the yield stress rises linearly 
with pressure and decreases monotonically with k s . The 
value of a drops by up to 50% as k s decreases from infin- 
ity to 210ecr~ 3 . Note that for k s = 210ecr~ 3 the rms dis- 
placement about lattice sites is near the usual Lindemann 
criterion for meltingj^j so much weaker values would be 
unphysical. To minimize computer time and the number 
of free parameters, we used rigid walls fc s — * oo in the 
following sections. 
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FIG. 6. Yield stress vs. pressure for the indicated values 
of the stiffness k B of the springs coupling wall atoms to lattice 
sites. The top panel shows results for — 90° (Fig. ^C) and 
the bottom panel shows results for 8 = 0° but d/d! = 12/13 
(Fig. |l|D). In each case, decreasing the stiffness decreases the 
yield stress. These results are from the search algorithm with 



1/2 coverage, o\ 
n=6. 



f = a, e w f = e, r, 



= 2 



rfnn = 1.2cr, and 
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E. Effect of wall orientation, coverage and chain 
length 

Calculated values of the static friction-fesibare surfaces 
vary dramatically with wall geometry.Bt3'E3 The yield 
stress equals that of a bulk crystal for identical aligned 
surfaces (6 — 0°) and is identically zero for incommen- 
surate orientations. Most experiments do not control 
the wall orientation, and yet measured friction coeffi- 
cients between nominally dry surfaces vary by roughly 
10% within ptiie-, same laboratory and by 20% between 
laboratories E3l£3 Another factor that is not controlled 
in most experiments is the density and chemistry of 
adsorbed molecules, although humidity is controlled in 
some cases. If adsorbed molecules are crucial in deter- 
mining the static friction, the friction they produce must 
be relatively insensitive to the wall orientation, sliding 
direction, coverage and chain length. 

Figure shows the yield stress as a function of pres- 
sure for the different walls shown in Figure ^. At low 
coverages (Fig. @(a)), the friction is nearly identical for 
all incommensurate walls. The friction for commensurate 
walls is roughly four times higher, but much less than the 
value for bare walls. At a coverage of 1/2 (Fig. 0(b)), the 
difference in the yield stresses for different incommensu- 
rate walls increases to roughly 20%. A similar variation 
with the direction of sliding (<fr = or 90°) is seen. Most 
of this variation is due to changes in To and all the curves 
are nearly parallel. Fit values of a, which dominates the 
high pressure friction, only change by about 5%. 

The greater variability at coverage of 1/2 arises be- 
cause this case corresponds to a dense monolayer between 
the surfaces. Molecules can not pick and choose which 
portions of the surface to cover. As a result the friction is 
higher for geometries B and D, where there are larger re- 
gions where the walls look commensurate and molecules 
can easily lock in registry with both surfaces. The fric- 
tion is smallest for geometry C, where there is almost no 
region where the two walls are in registry. 

Figure H shows the yield stress vs. pressure at different 
coverages for the extreme case of geometry G {9 = 90°). 
Note that the yield stress exhibits the usual linear pres- 
sure dependence at all coverages, and that most of the 
variation with coverage is in the intercept To. The value 
of a stays within 15% of the mean value, 0.030. At each 
pressure, the friction decreases as the coverage rises from 
1/8 to 1/2. During this increase in coverage, more and 
more of the sites with poor registry must be occupied. 
The friction goes down to a minimum value at a com- 
plete monolayer which corresponds to about 1/2 cover- 
age. If the difference between the sizes of the monomers 
and wall atoms was bigger, locking would be frustrated 
even in the areas with good registry. This might lead to 
even smaller variations with coverage than those shown 
here. 
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FIG. 7. Shearing direction and orientation dependence for 
(a) 1/8 coverage and (b) 1/2 coverage. The symbol shape 
indicates which of the geometries in Fig. [I] were used (see 
legend) . Open and closed symbols are for lateral force applied 
along the x (<f) — 0°) and y ((f) = 90°) directions, respectively. 
The ramp algorithm was used in (a) and the search algorithm 
was used in (b). The potential parameters were cr w f = a, 
e w f = e i r c = 2 1//6 <7, d nn = 1.2cr, and n=6. 



As coverage rises from 1/2 to 1 (about two monolay- 
ers) , the friction rises once more. The reason is that two 
layers of molecules have more internal degrees of freedom 
and can lock more easily into simultaneous registry with 
both walls. For this range of coverages yield continues 
to localize at the surfaces of the walls. Previous work on 
thicker films shows .that yield can move into the film as 
thickness increases.BEj In this limit r s approaches the 
bulk yield stress of the material making up the adsorbed 
layer, and eventually becomes independent of the cover- 
age. Note that the bulk yield stress of polymeric materi- 
als is also found, to follow the linear pressure dependence 
of Equation [l].El 

Figure ^ illustrates the interplay between the chain 
length and coverage dependence. At small coverages 
(Fig. H (a)), there is little dependence on chain length or 
wall orientation. For a complete monolayer (Fig. ^| (b)), 
chains with n — 3 and 6 show nearly the same behavior, 
while spherical molecules (n = 1) show a dramatically 
reduced friction and a. This suggests that the spheres 
are more constrained by their neighbors from locking into 
registry with the walls. This may reflect the fact that al- 
though all films contain the same number of monomers 
at a coverage of 1/2, the spheres are closer to their true 
equilibrium monolayer density than the chains. The rea- 
son is that the strong bonds within chain molecules lower 
the equilibrium spacing and thus increase the equilibrium 
density at a given pressure. The presence of two diffcr- 
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ent types of bonds in chain molecules may also improve 
their ability to rearrange to lock into both surfaces. An- 
other possibility is that chains produce a greater fric- 
tion because of the strong intramolecular bonds. If one 
monomer in the chain is locked to the bottom wall and 
another is in registry with the top wall, the intramolecu- 
lar bonds will lock the two walls together. This process 
has no analog for monomers and is particularly important 
as the coverage increases beyond a monolayer. 
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FIG. 8. Yield stress vs. pressure for the indicated coverages 
with system C (0 = 90°) where the largest variation was ob- 
served. Note that a coverage of 1/2 on each separated surface 
produces about a monolayer when the surfaces are placed in 
contact. Results are from the search algorithm with cr w f = cr, 



>V6, 



1.2a, and n—6. 



Figure |lC] shows the static friction as a function of 9 at 
a single pressure, P = 24eer~ 3 . This large pressure was 
chosen to maximize the contribution of a to the yield 
stress. Only 9 between and 30° are shown because 
other angles are related by symmetry. 

The commensurate case 9 — 0° shows a much larger 
friction than any of the incommensurate systems. Even 
a 2° rotation reduces the friction to a value that is com- 
parable to those at other values of 9. The same sharp 
transition was seen for coverages from 1/8 to 1. The data 
shown is for coverage of 1/2 where the variation of the 
yield stress at incommensurate 9 has already been shown 
to be largest. Even in this case the variation is less than 
±20%, which is comparable to experimental variations. 
The smaller variations of yield stress with 9 at a coverage 
of 1/8 are comparable to our statistical errorbars. 

The sliding direction 4> also has relatively little influ- 
ence on the yield stress between incommensurate sur- 
faces. Two orthogonal directions <f> = 0° and 90° are 
shown in Fig. |l^. The difference between them is most 
pronounced at small 9, where = 0° corresponds to 
pulling nearly parallel to the nearest-neighbor direction 
in each wall, and <fi = 90° to pulling along the perpendic- 
ular direction. 
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FIG. 9. Yield stress vs. pressure for the indicated values 
of chain length n at coverage of (a) 1/8 and (b) 1/2. Open 
(closed) symbols are for system C (B). The potential param- 
eters were a w f = a, e w f = e, r c = 2 1 / 6 er, and d nn = 1.2a. 
Errorbars are comparable to the system size and the ramp 
(search) algorithm was used for the top (bottom) panel. Ex- 
cept for spherical molecules at 1/2 coverage, the friction is 
insensitive to chain length. 

Figure [ll] shows the variation with <j> in more detail 
for coverage of 1/2 and 9 = 10.4°, where the variation in 
Fig. [l(] is large. The data are symmetric about 9/2 which 
is the direction of a 2-fold axis. This direction gives the 
maximum friction, and the yield stress drops to a plateau 
value on either side. The variation is within our statisti- 
cal errors at low pressures and is less than ±20% at the 
largest pressures studied. As in Figures 0(b) and |l^, slid- 
ing along the x axis leads to higher friction than sliding 
along y (equivalent to <p — and 30°, respectively). 

Note that when the yield stress is exceeded the top 
wall does not always move exactly in the direction of 
the applied force. This effect is most pronounced for the 
commensurate case, but is also seen at small 9. 



F. Effect of Interaction Potentials 

Experimental results for different sliding surfaces and 
different surface layers show much larger variations than 
the fluctuations for a given material system. Such 
changes in material would correspond to changes in the 
parameters in our interaction potentials such as the 
strength of the interaction between the wall and fluid e w f , 
the characteristic length of this interaction <r w f , the range 
of the interaction r c , and the nearest-neighbor spacing 
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d nn in the walls. Other structural properties and the 
elastic moduli and hardness of the solids may also be 
important, but are not easily included within our model. 




FIG. 10. Yield stress as a function of orientation angle 
8 at P = 2Aea~ 3 . The commensurate case corresponds to 
6 — and the smallest non-zero angle shown here is 2.02°. 
These results are from the search algorithm with 1/2 cover- 
age, cr w f = a, e w f = e, r c = 2 1 ^ 6 a, d nn = 1.2a, and n=6. The 
simulation cell was 57.6a by 49.9<r. 




FIG. 11. The yield stress vs. the angle (f> between the ap- 
plied stress and the x axis for P = 12ea~ 3 (circles) and 32e<7 -3 
(squares). These results are from the search algorithm with 
1/2 coverage, cr w f = a, e w f = e, r c = 2 1 ^ 6 a, d nn = 1.2cr and 
n = 6. The simulation cell was 57.6<r by 49.88a. Errorbars 
are only shown when they are larger than the symbol size. 



The results shown so far were all for purely repulsive 
potentials, r c = 2 1 / 6 cr. Increasing the cutoff to r c = 2.2a 
includes most of the attractive tail in the Lennard-Jones 
potential. This attractive tail leads to an adhesive force 
pulling the walls together in addition to the external pres- 
sure P. As shown in Fig. 12(a) and (c), the yield stress 
for r c = 2.2a is parallel to results for r c = 2 1 / 6 cr and 



could be obtained by shifting the zero of pressure to re- 
flect the adhesive stress. The required shift increases 
with increasing coverage because the density of adsorbed 
molecules that contribute to the adhesive interaction in- 
creases. 

When the above results are fit to Eq. [j], the value of 
a is independent of r c , and the value of To increases with 
the amount of adhesion. For purely repulsive interactions 
To is negative because a finite pressure is n eeded to lock 
the adsorbed layers in a glassy state (Sec. Ill Op . When 
r c is increased to 2.2a, a glass forms at zero pressure, 
and the value of To is positive. Many authors refer to To 
as the adhesive contribution to friction, and it is obvi- 
ous that to I a represents the effective adhesive pressure 
that must be added to the external pressure to get strict 
proportionality between the friction and total load. 

These simulations did not include the direct interac- 
tions between atoms from different walls. As noted in 



Section [I A , the wall spacing is large enough that the 
direct interaction vanishes for Even for 

r c = 2.2er, wall interactions have little effect. For exam- 
ple the extra adhesive pressure due to wall/monomer in- 
teractions is To/a ~ 5eer~ 3 for 1/2 coverage (Fig. |l2|(a)). 
Direct calculation of the extra adhesive pressure due to 
interactions between atoms on different walls yields a 
value that rises from 0.4 to 1.3ec~ 3 over the observed 
range of wall separations. Correcting the external pres- 
sure to account for this interaction would change a by 
less than 3%. 

Panels (b) and (c) of Fig. |l2| show that doubling 
the strength of the potential between the surfaces and 
monomers has almost no effect on the yield stress. The 
main factor that changes a is the ratio between the char- 
acteristic separation between wall and fluid atoms, cr w f, 
and the nearest-neighbor spacing between wall atoms, 
d nn . The value of a rises monotonically as a w f/d nn de- 
creases. For the cases shown, a changes by an order 
of magnitude. Smaller values of a are obtained by in- 
creasing er w f /dim further, but it is difficult to increase 
a for this geometry. The reason is that the adsorbed 
molecules begin to penetrate into the wall when a w f/d nn 
is too small. Other crystallographic orientations of the 
surface, such as (100), may produce larger values of a. 
Flat buljdisordered surfaces produce values of a as high 
as 0.3.ffl 

The strong effect of a w { /d nn and weak effect of e w f / e 
suggest a geometric origin of a that is analogous to 
early geometric theories for AmptLtons' lawsEJ and to Is- 
raelachvili's cobblestone modelJla'E 3 ] At pressures where 
the adsorbed layer acts like a glass, the repulsive part of 
the Lennard-Jones interaction dominates. The repulsive 
force between monomers and wall atoms at a separation r 
is given by F rep = 48(e w f/<7 w f)(cr w f/r) 13 . The mean value 
of F Tep must scale with the pressure. Hence, increasing 
P by a factor of 2 only decreases the typical separation 
between monomers and wall atoms by a factor of 2 1 / 13 or 
5%. Increasing e w f by a factor of two leads to the same 
small decrease in separation. In contrast, the separation 
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changes by a factor of 2 12 / 13 w 1.9 with a factor of 2 
change in <r w f. Thus monomers and wall atoms can be 
thought of as hard spheres whose diameter D is nearly 
independent of pressure and e w f but directly proportional 

to <7 w f. 

Figure [l3] illustrates the effect of changing the ratio be- 
tween the hard sphere diameter and the nearest-neighbor 
spacing d nn . The hard-sphere repulsion leads to a sur- 
face of closest approach between a monomer and the wall. 
In order for the monomer to slide relative to the wall it 
must be lifted up the ramp formed by this surface. The 
lateral force required to move monomers up the ramp is 
just the normal load times the slope of the surface. Thus 
a will be proportional to some average of the slope over 
all monomers. As D/d nn decreases, monomers penetrate 
more deeply between the surfaces and both the slope and 
a increase. This is precisely the trend seen in Fig. [l^. 
Thus the hard-sphere model explains the linear depen- 
dence of t on pressure as well as the insensitivity to e w f 
and the variation with er w f and d nn . 
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FIG. 12. Yield stress vs. pressure for various choices of the 
potential parameters in system C with n — 6. The coverage 
is 1/2 in (a) and (b) and 1/8 in (c). Unless noted the default 
values of <r w f = a, e w f = e, d nn = 1.2<r, and r c = 2 1 / 6 <r are 
used. Each of these parameters was changed in turn and the 
corresponding results are labeled by the changed value. The 
search algorithm was used for the data in (a) and (b), and 
the ramp algorithm was used in (c). 
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FIG. 13. Sketch of surface of closest approach for 
D/dnn = 0.7 (solid line) and 1.3 (dashed line). The lateral 
force required to lift monomers up the ramps that these sur- 
faces represent is given by the slope times the normal force. 



IV. SUMMARY AND CONCLUSIONS 

The results presented in this paper explore the effect of 
a specific type of third-body small adsorbed molecules, 
on the friction between the bounding solids. The effect 
of temperature, coverage, crystalline alignment, lattice 
constant, system size, molecular size and interaction po- 
tential parameters were all examined. 

The first important result is that static friction is ob- 
served for all systems - even when the static friction be- 
tween bare walls is identically zero. The shear stress 
needed to initiate sliding approaches a constant value in 
the thermodynamic limit (Fig. 0) and follows the linear 
pressure dependence (Equation Th) seen in many exper- 
imental systems. Moreover, the shear stress is insensi- 
tive to parameters that are not controlled in most ex- 
periments, including crystalline alignment, coverage and 
sliding direction. Even under conditions where the vari- 
ations with these parameters are largest, the changes are 

jariations between 
''^ This is an impor- 



comparable to typical experimen 
different samples and laboratories 
tant test of any molecular scale model for friction. 

The value of tq in Equation [l] is increased by attractive 
interactions between the monomers and walls or between 
the two walls. The adhesive interaction was increased 
both by increasing coverage (Fig. ||) and by increasing 
r c to increase the contribution from the attractive tail 
in the Lennard- Jones potential (Fig. |l2|). Increasing the 
temperature (Fig. |^) increases the entropic repulsion and 
lowers To. 

For each system there is a pressure Pt — —r^ja be- 
low which the adsorbed layer acts like a liquid lubricant 
(Sec. Ill C ) . When the attractive tail of the Lennard- 
Jones potential is included (large r c ), the value of P t 
for submonolayer films is negative up to quite high tem- 
peratures. This suggests that thicker films are needed 
to achieve lubricated behavior in typical experimental 
systems. However, it would be interesting to explore 
whether the static friction between flat surfaces can be 
eliminated by the proper choice of molecule and the use 
of elevated temperatures. 

All of the observed trends in a can be understood in 
terms of a simple geometrical model where monomers act 
like hard spheres. Within this model, a only depends on 
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ile atoms 
rearrange 



the hard sphere diameter and the wall geometry. Param- 
eters such as pressure, temperature, wall/monomer cou- 
pling energy, and r c have almost no effect on the hard 
sphere diameter and thus a. In contrast, changes in cr w f 
or d nn produce direct changes in the geometry and change 
a dramatically (Fig. |l2| ). Smaller shifts in a are pro- 
duced by changes in coverage that alter the ability of the 
hard spheres to pack into a configuration that locks into 
both walls (Fig. |). 

Simple model potentials were used for the adsorbed 
molecules and surfaces in this study, in order to allow a 
wide range of parameters to be investigated. More re- 
alistic models will have more complex interactions and 
multiple atomic species. This will complicate the geom- 
etry, but one expects that the basic idea of an effective 
surface of closest approach defined by hard-sphere repul- 
sions will still be valid at the pressures of interest. This 
naturally produces static friction and a shear stress that 
obeys Equation [l]. Indeed the same picture can be ap- 
plied to larger particles between bearing surfaces such as 
dust, sand or wear debris. 

While our geometric argument for a bears some 
similarity to-.earlier models, based on interlocking of 
macroscopicE3 or atomicE^rtJ asperities, there is an im- 
portant difference. These previous models envisioned as- 
perities that were rigidly attached to their respective sur- 
faces and would translate with them, 
considered here and in previous work 
to create interlocking configurations for any relative po- 
sition of the surfaces. Not only does this explain why 
static friction can occur between surfaces with uncorre- 
cted asperities, it also makes it possible to connect static 
and kinetic friction. _ 

As pointed out by LeslieEll in 1804, interlocking of rigid 
asperities can explain static friction, but predicts vanish- 
ing kinetic friction. The reason is that the energy needed 
to move up ramps is recovered as the asperities move back 
down and the total work required to slide the surfaces 
vanishes. The situation is very different when the inter- 
locking occurs dynamically through motion of adsorbed 
molecules. As shown in recent workJ23'C3 at any given 
instant almost all of the mobile molecules are trapped 
in local free energy minima. Since lateral motion of the 
walls increases their energy within this minimum they 
contribute to the friction. The kinetic and static friction 
are closely related because the distribution of molecules 
in free energy minima is nearly the same. The differ- 
ence is that when the walls are moving at fixed velocity 
the minimum holding each molecule eventually becomes 
unstable. It then pops forward rapidly to the nearest 
metastable state. Much of its energy is dissipated in this 
pop because molecules pop forward at different locations 
at different times in an incoherent manner and the po- 
tential they move to varies in a random way. It is not 
possible for them to move up and down over a rigid po- 
tential energy surface as in the models of fixed asperities 
considered by Leslie and earlier researchers. 
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